Polymer translocation through a nanopore - a showcase of anomalous diffusion. 
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The translocation dynamics of a polymer chain through a nanopore in the absence of an external 
driving force is analyzed by means of scaling arguments, fractional calculus, and computer simula- 
tions. The problem at hand is mapped on a one dimensional anomalous diffusion process in terms 
of reaction coordinate s (i.e. the translocated number of segments at time i) and shown to be 
governed by an universal exponent a = 2j{2v + 2 — 71) where v is the Flory exponent, and 71 is 
the surface exponent. Remarkably, it turns out that the value of a is nearly the same in two- and 
three-dimensions. The process is described by a fractional diffusion equation which is solved exactly 
in the interval < s < N with appropriate boundary and initial conditions. The solution gives the 
probability distribution of translocation times as well as the variation with time of the statistical 
moments: (s(t)}, and (s 2 (t)) — (s(t)) 2 which provide full description of the diffusion process. The 
comparison of the analytic results with data derived from extensive Monte Carlo (MC) simula- 
tions reveals very good agreement and proves that the diffusion dynamics of unbiased translocation 
through a nanopore is anomalous in its nature. 

PACS numbers: 82.35.Lr, 87.15.Aa 



The dynamics of polymer translocation through a pore 
has recently received a lot of attention and appears highly 
relevant in both chemical and biological processes [lj] . The 
theoretical cosideration is usually based on the assump- 
tion 0, 0, IH that the problem can be mapped onto a one 
dimensional diffusion process. The so called transloca- 
tion coordinate (i.e., reaction coordinate s) is considered 
as the only relevant dynamic variable. The whole poly- 
mer chain of length A is assumed as being in equilibrium 
with a corresponding free energy J-(s) of an entropic na- 
ture. The ld-dynamics of the translocation coordinate 
follows then the conventional Brownian motion and the 
one dimensional Smoluchowsky equation [B| can be used 
with the free energy F(s) playing the role of an external 
potential. In the absence of external driving force (un- 
biased translocation), the corresponding average first - 
passage time follows the law r(A) oc a 2 N 2 /D, where a 
is the length of a polymer Kuhn segment and D stands 
for the proper diffusion coefficient. The question about 
the choice of the proper diffusion coefficient, D, and the 
nature of the diffusion process, is controversial. Some 
authors @, 0] adopt D oc A -1 , as for Rouse diffusion, 
which yields t oc A 3 as for polymer reptationQ albeit 
the short pore constraint is less severe than that of a tube 
of length A. In ref. [J] it is assumed that D is not the dif- 
fusion coefficient of the whole chain but rather that of the 
monomer just passing through the pore. The unbiased 
translocation time is then predicted to go as r oc A 2 . The 
latter assumption has been questioned 0| too. Indeed, 
on the one hand, the mean translocation time scalesQ as 
r ~ A 2 but on the other hand the characteristic Rouse 
time (i.e. the time it takes for a free polymer to diffuse 
a distance of the order of its gyration radius) scales as 



trousc oc N 2 " +1 where the Flory exponent v = 0.588 at 
d = 3, and v = 0.75 at d — 2 Q. Thus tr ousc 3> r against 
common sense given that the unimpeded motion should 
be in any way faster than that of a constrained chain! 
Moreover, the equilibration of the chain is questionable 
when the expression for J-(s) is to be used. The char- 
acteristic equilibration time scales again as r eq oc N 2v+1 
and is thus always larger than the translocation time, i.e. 
i~eq 3> t. Again the internal consistency of the whole a 
proach is in doubt. It was found by MC-simulation [~ 
that r oc A 2 - 5 for translocations in d = 2. This indicates 
that the translocation time scales in the same manner 
as the Rouse time albeit with a larger prefactor that de- 
pends on the size of the nanopore. The authors argued 
that this finding bears witness of the failure of the Brow- 
nian nature of the translocation dynamics and suggested 
instead that anomalous diffusion dynamics Q should be 
more adequate. The r oc A 2 ' 5 scaling law has been cor- 
roborated by a further MC - stu dy [1C| as well as by a 
MC-simulations on a 3d - lattice [llj . and it was shown 
that r oc jV 2 - 46±0 03 . The time variation of the second 
statistical moment, (s 2 ) — (s) 2 oc t a , clearly indicates 
an anomalous nature llj since the measured exponent 
a = 0.81 ±0.01 while r oc N 2 l a . Still missing is a proper 
theoretical analysis which could explain the physical ori- 
gin of the anomalous dynamics, and make it possible to 
solve the appropriate fractional diffusion equation (or, 
in case of a biased translocation, the fractional Fokker - 
Planck - Smoluchowski equation) 0, [l2| governing this 
dynamics. 

In this Letter we suggest a unique physical picture 
which justifies the mapping of the id problem on a Id re- 
action cooordinate, s, and we show that the latter obeys 
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FIG. 1: Schematic representation of a chain fold of length s 
moving through a nanopore. The transition rate is slowed 
down by an entropic barrier: a) initially the fold is on the cis- 
side of the wall, b) the fold entropy decreases during threading 
because of the fold fragmentation. 



anomalous diffusion dynamics, described by a fractional 
diffusion equation. We solve this equation exactly, sub- 
ject to the proper boundary conditions, and find a per- 
fect agreement with our scaling prediction. Eventually, 
we demonstrate that the results of our 3d off-lattice MC- 
simulations are in accord with our analytical findings. 

Mapping onto Id dynamics - As indicated above, the 
assumption that the whole polymer chain is in equilib- 
rium and the diffusion is governed by conventional Brow- 
nian dynamics leads to contradictions. Instead, we as- 
sume now that only a part of the whole chain may equili- 
brate between two successive threadings. This part of the 
chain which adjoins the membrane on the cis- or trans- 
sides will be denoted as fold and we assume that it is 
much shorter than the whole chain length N but is still 
long enough so that one can use the principles of statis- 
tical physics. We also assume that the excluded volume 
interaction of a fold with the rest of the chain is relatively 
weak so that it could be treated as a subsystem with a 
well defined free energy. This latter assumption is based 
on the observation that the chain on either of the two 
sides may be seen as a polymeric "mushroom" whereby 
the monomer density close to the membrane (or wall) is 
much smaller than the density inside a single coil (see 
Fig.4 and eq.(II.4) in ref.[I3]). Thus one can claim that 
there is a depletion area near the membrane [l4T | 

Figure[T]illustrates how a fold squeezes from the cis- to 
the trans - side through a short nanopore (of length « a) 
which is slightly wider than the chain itself. It is self - 
evident that, in the absence of the external force, with the 
equal probability folds from the trans- side could go to the 
cis - side. If the irons-part of the fold in FigQ]has length 
n then the corresponding free energy function F t (n)/T — 
—n In k — (71 — 1) In n, where n is the connective constant 
and 71 is the surface entropic exponent [H| . For the cis- 
part of the fold one has F c (n) /T = — (s — n) In k — (71 — 
l)ln(s — n) so that the total free energy is T(n)jT = 
— s InK — (71 — 1) ln[n(s — n)]. One can, therefore, ascribe 
to the fold cis-trans - transition a pretty broad barrier 
given by F{n). The corresponding activation energy of 
the fold can be calculated as AE(s) = F{s/2) - = 
(l- 7 i)Tln«. 



How can we estimate the characteristic time for the 
fold cis - trans - translocation? In the absence of a sep- 
arating membrane this would be the pure Rouse time 
£r. oc s 2u+1 0. The membrane with a nanopore imposes 
an additional entropic activation barrier AE(s) which 
slows down the transition rate. The characteristic time, 
therefore, scales as t(s) = i R (s) exp[A_E(s)] oc s 2i/+2 ~ 71 . 
This makes it possible to estimate the mean-squared dis- 
placement of the s - coordinate: 



(s 2 ) cxt 2 ^ 2v+2 -^\ 



(1) 



Hence the mapping onto the s - coordinate leads to an 
anomalous diffusion law, (s 2 ) oc t a , where a = 2f{2v + 
2 — 71). Taking into account the most accurate values of 
the exponents for d — 3 , v — 0.588 and 71 = 0.680 [lj|, 
we obtain a — 0.801. In turn, the average translocation 
time r oc N 2 / a oc JV 2 496 . Remarkably, in 2d where V2d — 
0.75 and 71 w O.9450, one finds a w 0.783 (i.e., a 
is almost unchanged)! This explains why the measured 
exponents in both 2d @ and in 3d 11 1 are so close. The 
derivation of a is our central scaling prediction. It also 
agrees well (see below) with our own MC - data on the 
translocation exponent. 

Fractional diffusion equation - We now turn to the 
fractional diffusion equation (FDE) which furnishes a 
natural framework for the study of anomalous diffusion 

[l2| . Here we make use of this method in a systematic 
way. Our FDE reads 



^W(s,t) = Dl-*K a -^W(s,t) 



(2) 



where W(s, t) is the probability distribution function 
(PDF) for having a segment s at time t in the 
pore, and the fractional Riemann - Liouville opera- 
tor vDl- a W{s,t) = (l/T(a))(d/dt) f Q dt'W{s,t')/{t - 
t') 1 ^ a . In Eq. ^ T(a) is the Gamma-function, and 
K a is the so called generalized diffusion constant. This 
constant could be defined as K a = T(l + a)l 2 /(2r^) in 
terms of the fold length / and the waiting time scale r w 
(see Chapter 3.4 in [9(). It should be mentioned that 
the constant K a is the only adjustable parameter of our 
theory which will be fixed below through the comparison 
with our MC - data. 

Recently the method of generalized Langevin equation 
(GLE) has been used to describe anomalous conforma- 
tional dynamics within single molecule proteins [H}. In 
contrast to the FDE - approach which deals with the 
total distribution function at particular boundary con- 
ditions (see below), the GLE - method treats only the 
first two moments (or time-correlation functions, mem- 
ory kernel, etc.). To the best of our knowledge, at the 
present time it is not clear how one can derive in a closed 
form a non-Markovian Fokker-Planck equation for the 
distribution function [l9| starting from the GLE. On the 
other hand, the translocation time distribution function 
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(see below) is an entity of great importance because it 
could be directly measured in experiment Therefore, 
we prefer to use FDE - approach for the translocation 
problem. 

Consider the boundary value problem for FDE in the 
interval < s < N. This problem has been dis- 
cussed before in the context of even more general frac- 
tional Fokker-Planck equation (FFPE) 0. The bound- 
ary conditions correspond to the reflecting- adsorbing 
case, i.e., (d/ds)W(s,t)\ s = = and W{s = N,t) = 

0. The initial distribution is concentrated in sq-, i.e., 
W(s,t — 0) = 8(s — s ). The full solution can 
be represented as a sum over all eigenfunctions (/? n (s), 

1. e., W{s,t) = X^o^M^nfs) where ip n (s) obey the 
equations K a (d 2 /ds 2 )ip n (s) + X n , a ip n (s) = 0, and the 
eigenvalues X n , a can be readily found from the fore- 
going boundary conditions, as a result \ n , a = (2n + 
1) 2 tt 2 A q /(4A 2 ). The temporal part T n (t) obeys the 
equation (d/dt)T n (t) = -\ n , a vD]~ a T n (t). The so- 
lution of this equation is given by T n (t) — T n {t 
G)E a (-X n>a t a ) where the Mittag - Leffler func- 
tion E a {x) is defined by the series expansion E a (x) = 
J2^Lo x k /T(l + ak). At a = 1 it turns back into a stan- 
dard exponential function (normal diffusion). Thus we 
arrive at the complete solution of eq. ( [2]) : 
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First - passage time distribution - The distribution of 
translocation times (which could, in principle, be mea- 
sured in experiment) is nothing but the first - passage 
time distribution (FPTD) Q(so,t), where sq stands for 
the initial value of the s - coordinate. Knowing the prob- 
ability distribution function W(s,t), we can calculate 
FPTD Q(so,t). The relation between both functions is 
given as Q(s ,t) = -(d/dt) J* W{s,t)ds @|. This yields 
the FPTD as follows 
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a laboratory experiment there always exists some upper 
time limit t*. Taking this into account, one can show that 
an "experimental" first passage time scales as r ~ N 2 / a 
2l| which we observe in our MC-simulation. 

Statistical moments (s) and (s 2 ) vs. time - The sub- 
diffusive behavior of the second moment (s 2 ) — (s) 2 cx t a 
is a hallmark of anomalous diffusion. Starting from eq. 
([3]) we can immediately calculate them. The calculation 
of the first moment (s) 
yields 



sW(s,t)ds/ J" W(s,t)ds 



(s) (t) =i _2E„= 
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(5) 



Since E a (t = 0) = 1, the initial value (s)(t = 0) = 
(we put So = 0) as it should be. In the opposite 
limit, t — > oo, we can use the asymptotic behavior 
E a [— \ n a t a ] — l/r(l — a)X n a t a as well as the sums val- 
= ues £io l/(2n + l) 4 = ^ 4 /96 and E^("l)V(2n + 
l) 3 = 7r 3 /32 in the nominator and denominator respec- 
tively. After that (s)(t — > oo) = A/3, i.e. the function 
goes to a plateau. 

The result for the second moment, (s 2 ) = 

L s 2 W(s, t)ds/ W(s, t)ds, can be cast in the follow- 
ing form 
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(6) 



(4) 



where the generalized Mittag - Leffler function E a ^ a (x) = 
Efc=o xk /F{a + ka). The long time limits of Mittag - 
Leffler functions in eqs. ((3]) and |4]) follow an inverse 
power law behavior, E a (— X n . a t a ) oc l/r(l — a)X niCc t a 
and E a _ a (— X n ^ a t a ) oc a/T(l — a)X 2 l a t 2a . By making use 
of this in eq. (gj), the long time tail of the FPTD then 
reads Q(t) oc aA 2 /2r(l - a)K a t 1+a . This behavior is 
checked below in our MC - investigation. It can be seen 
that the mean first passage time, defined simply as r = 
J °° tQ(t)dt, does not exist 21, 22|. On the other hand, in 



Again it can be readily shown that (s 2 )(0) — (s) 2 (0) = 0. 
In the long time limit in the same way as before and tak- 
ing into account that ErT=o(- 1 )'V(2" + I) 5 = 5tt 5 /i536 
we will find (s 2 ) (t -> oo) - (s) 2 (*-►<») = A 2 /9 

Monte Carlo data vs. theory - We have carried out 
extensive MC - simulations in order to check the main 
predictions of the foregoing analytical theory. We use 
a dynamic bead-spring model which has been described 
before [24J , therefore we only mention the salient features 
here. Each chain contains A effective monomers (beads), 
connected by anharmonic FENE (finitely extensible non- 
linear elastic) springs, and the nonbonded segments in- 
teract by a Morse potential. An elementary MC move 
is performed by picking an effective monomer at random 
and trying to displace it from its position to a new one 
chosen at random. These trial moves are accepted as 
new configurations if they pass the standard Metropolis 
acceptance test. It is well established that such a MC 
algorithm, based on local moves, realizes Rouse model 
dynamics for the polymer chain. In the course of the 
simulation we perform successive run for chain lengths 
N = 16, 32, 64, 128, 256 whereby a run starts with a con- 
figuration with only few segmens on the trans-side. Each 
run is stopped, once the entire chain moves to the trans- 
side. Complete retracting of the chain back to the cis-side 
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N=16 




FIG. 2: The translocation time distribution function Q(t): a) 
Scaling plot of the theoretical predictions calculated from Eq. 
l|4|l for different chain lengths N. Dashed line denotes the 
long time asymptotic tail with slope —1.8. The inset shows 
Q(t) for N = 256 in normal coordinates, b) The FPDT Q(t) 
from the MC-simulation for N = 64, 128. The inset shows the 
expected (r) vs. chain length N dependence and the straight 
line is a best fit with slope ~ 2.52 ± 0.04. 



is prohibited. During each run we record the transloca- 
tion time t, and the translocation coordinate s(t). Then 
we average all data over typically 10 4 runs. In principle, 
the pore may apply a drag force on the threading chain 
due to a chemical potential gradient, however, in the 
present work we consider only unbiased diffusion. In Fig- 
ure [2] (a) we show a master plot of the translocation time 
distribution Q(t) derived from Eq. Q for different chain 
lengths N = 16, 32, 64, 128, 256. For the calculation of 
data we have used Mathematica with a special package for 
computation of Mittag-Leffler functions [5^1 . Evidently, 
all curves collapse on a single one when time is scaled as 
t oc N 2 l a with the predicted a = 0.8. The long-time tail 
for this value of a should exhibit a slope of —1.8. The 
inset in Figure [UJa) reveals the long tail of Q(t) for large 
times. A comparison with Figure^b) demonstrates good 
agreement with the simulation data despite some scatter 
in the FPDF even after averaging over 10000 runs. As 
shown in the inset, the mean translocation time scales as 
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FIG. 3: Variation of the first and second moments of the PDF 
W(s,t) with time for chain lengths N = 16,32,64, 128,256. 
a) Log-log plot of the first moment (s) vs. time t from a MC 
simulation (big symbols) and from Eq. © (small symbols). 
The dashed line denotes t a ^ 2 with a slope of 0.4. b) The same 
as in a) but for (s 2 ) — (s) 2 . Analytical data is obtained from 
Eq. ©. The dashed line has a slope of 0.8. 



(t) oc N 2 - 5 in good agreement with the predicted a = 0.8. 
An inspection of Figure [3] where the time variation of the 
PDF W(s, t) moments are compared demonstrates again 
that data from the numeric experiment and the analytic 
theory agree well within the limits of statistical accuracy 
(which is worse for N — 256). Not surprisingly, the time 
scale of the MC results does not coincide with that of Eqs. 
([5]), ([6]) since in the latter we have set K a , which fixes the 
time scale, equal to unity. Closer examination of Figure 
[3] shows that the resetting of the generalized diffusion co- 
efficient as K a ~ (80) 8 ~ 33.3 enables to superimpose 
the results of theoretical calculation and MC - data. 

In summary we have shown unambiguously that the 
translocation dynamics of a polymer chain threading 
through a nanopore is anomalous in its nature. We 
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have succeeded to calculate the anomalous exponent 
a = 2/(2v + 2 — 71) from simple scaling arguments and 
embedded it in the fractional diffusion formalism. Wc 
derived exact analytic expressions for the translocation 
time probability distribution as well as for the moments 
of the translocation coordinates which are shown to agree 
well with our MC simulation data. The present treat- 
ment can be readily generalized to account for a drag 
force on the chain and results for this case will be re- 
ported in a separate publication. 
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